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We investigate the Dirac time-dependent variational method for a system of non-ideal 
Bosons interacting through an arbitrary two body potential. The method produces a set 
of non-linear time dependent equations for the variational parameters. In particular we 
have considered small oscillations about equilibrium. We obtain generalized RPA equations 
that can be understood as interacting quasi-bosons, usually mentioned in the literature as 
having a gap. The result of this interaction provides us with scattering properties of these 
quasi-bosons including possible bound-states, which can include zero modes. In fact the zero 
mode bound state can be interpreted as a new quasi-boson with a gapless dispersion relation. 
Utilizing these results we discuss a straightforward scheme for introducing temperature. 
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1 Introduction 



The objective of this paper is to exhibit the most general way of obtaining time dependent 
equations of motion in the Gaussian approximation This will lead to the so called 
generalized RPA, when one examines infinitesimal oscillations about equilibrium. The static 
solution in the uniform case can be obtained using several other methods - [[| leading to a 
gap in the quasi-boson energy. We show here that the time dependent RPA equations lead to 
a gapless mode. In fact this must happen because particle number conservation symmetry 
was broken in the static solution, so the zero gap is exactly the associated Goldstone mode. 
This can be seen as an alternative to the functional derivative [|J method in the Girardeau- 
Arnowitt || approximation. We stress that to obtain this result the equations of motion 
must arise from a consistent variational scheme. An important element of this scheme is the 
fact that even in the uniform case the variational quantities have to depend on two momenta. 
It is actually this feature which leads to a gapless result. Since this solves the problem of 
the gap, truncations normally used to avoid it || are not necessary. 

We would like to point out that due to the likelihood of increased densities in future 
experiments on Bose-Einstein condensation in atomic traps - |§ our time dependent 
equations may be useful to study non-equilibrium evolution. The RPA equations also give us 
the possibility of investigating the scattering of two quasi-bosons in the medium. Information 



about scattering can be used in calculating the free energy at finite temperature [|10] which 
takes into account the gapless dispersion relation in contrast to the so called Temperature 
Hartree Fock Bogoliubov method [13| . 

The structure of this paper is as follows. Section 2 reviews the time-dependent variational 
principle and the canonical nature of the equations of motion arising from it. In section 3 
we discuss the static solution. Section 4 obtains the small oscillation approximation (RPA) 
for the equations of motion. The Goldstone mode is discussed in Section 5. In section 6 
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we show the connection with the Bogoliubov transformation and calculate the one and two 
quasi-boson wave functions. In section 7 we show how the two quasi-boson scattering affects 
the thermodynamics using well know methods os statistical mechanics. 

2 General Formalism 

In this section we shall review [] Time- dependent Variational Principle and show how 



it can be implemented in the non-relativistic many-body case. First we define an effective 
action functional for the time-dependent quantum system 



S = J L(t)dt = J dt(V,t\(id t - H)\V,t), 



X 



where t) is the quantum state of the system and H is the Hamiltonian of the theory. For 
a system of non-relativistic interacting Bosons we have [we use the notation : J x = / d 3 x] 

H=[ ^(x)t/ l (x,y)^(y) + i / ^(y)^(x)V(x - y)$(x)^(y), (2) 

Jx,y I Jx,y 

where the one body Hamiltonian ft,(x, y) may include a one body external potential. The 
operators ip and $ can be written in the form 

$(x) = -= U(x) + ztt(x)1 (3) 



^(x) f = —= U(x) - ztt(x)1 (4) 



y/2 
1 

V2 

where 0(x) is the field operator and 7r(x) the canonical field momemtum 

We can obtain the time dependent Schrodinger equation by requiring that S is stationary, 
supplemented by appropriate boundary conditions, under the most general variation of t). 
The variational scheme is implemented by chosing a trial wave functional describing the 
system. Working in the functional Shrodinger picture we replace the abstract state |\&, t) by 
a wave functional of the field 0'( x ) 

l*,*)-*^]. (5) 



1 This material is substantially the same as some sections of JT3] and it is included here for completeness 
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The action of the operators 0(x) and the canonical momentum 7r(x) are realized respectively 
by 



0(x)|tt,t) - 0'(x)*[0',t] 
7T(X)|*,*) - * 



(6) 
(7) 



The mean value of any operator is calculated by the functional integral 

{%t\o\%t) = J \v^*w ,t\mw \t\, (8) 

where \& is normalized to unity. The Gaussian approximation consists in taking a Gaussian 
trial wave functional in its most general parametrization 

'G-\x,y,t) 



V[<j>',t] = Nexpl - f 5<f)'(x,t) 

+ z£7T(x,t)^ , (x,t)|, 



i£(x,y,t) 



(9) 



with S<f>'{x,t) = 0'( x ) — 0( x )^)- Due to the fact that the Hamiltonian commutes with the 
number of particles 

[H,N] = (10) 



we can define a more general trial functional 



\*',t) = e- i * 9 ®\*,t), 



(11) 



where 6{t) is another variational parameter introduced because of this continuous symmetry. 
Thus our variational parameters are 0(x, t), 7r(x, t), 9(t),G(x, y, t) and S(x, y, t) ,with G and 
E being real symmetric matrices. These quantities are related to the following mean-values: 



(*',i|0(x)|*',*> 
(tt',f|7r(x)|tt', t ) 
(*',*|0(x)0(y)|*',t) 



0(x,t) 
7r(x,t) 

G(x,y,t) + 0(x,t)0(y,t) 



(12) 
(13) 
(14) 



(*',t|7r(x,t)7r(y,t)|vl/',t) = G 1( *' y>t) +4/ S(x, w, f)G(w, z, f)E(z, y, t) 

4 Jw,z 

+7r(x,t)7r(y,t) (15) 
i« t|0(x, t)7r(y, t)|tf, t> = JJ*-ll + i jf [E(x, z, t)G(z, y, t) + G(y, z, t)E(z, x, t)} 

+20(x,t)vr(y,t) (16) 
= / E(x,y,t)G(y,x,t) + [ n(x,t)j>(x,t) 

+J\f9(t) + total time derivatives. (17) 

We may ignore the total time derivatives because they do not contribute to the equations of 
motion. If now we write the action we will get 

S = l dt {i 7r(X ' ^ (X ' t)+ iy S(X ' y ' t)6(y ' X ' ^ + N ® it) ~ U ) ' (18) 

where 

H = (ty',t\H\ty',t) (19) 

and 

7V= (V',t\N\V',t). (20) 

From (fl8|) we see that (J\f, 9), (tt, <fi) and (E, G) are canonical pairs. Because of the symmetry 
7i has no dependence on 9 and it follows that N = and 9(t) = constant = /i . We can now 
write the remaining Hamilton equations, 

5(H - MQ 

0(x ' t} = ^) ' (21) 

^ = - j ^mp (22) 

s(x , y , () = JJ^O. (24) 



For convenience we introduce 



, ? . 6(x, t) + Z7r(x, 
^(x, t) = (V>(x)> = n 1 ; , (25) 
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so that Eqs.@ID-(f22D become 



5(H- nM) 



It is well known that a Gaussian functional leads to the mean field factorization i.e 
H-\iN= / x , y {[Mx,y)-/i5(x-y)]p(x,y ) t) + ^(x-y)|^(x,t)| 2 |^(y,t)r} 



(26) 



and 



+\ L, y V(x - y) [R(y, x, t)R(x, y, t) + i?(x, x, t)i?(y, y, f) + L>*(x, y, t)D(x, y, t)] 



+ / Xiy V(x-y) 



^*(x,^(y,t)i?(x,y,t) + ^V*(y, t)V(x, t)i?(y, x, t) + |</>(x, t)| 2 i?(y, y, t) 



/ x , y F (x - y) [V>(x, t)^(y, t)£>*(x, y, t) + ^*(x, t)^*(y, t)D(x, y, t)] 



(27) 



p(x,y,t) = (^(x)^(y))=f( X ,^(y,t) + J?(xj,i) 
A(x,y,t) = -(V(x)^(y)) =--0(x,t)V(y,t)+D(x,y,t) 



(28) 
(29) 



with 



#(x,y,t) 



^(x,y,t) 



G^fry,*) 



+ G(x,y,t)-«J(x-y) 



2 / E(x,w,t)G(w,z,t)E(z,y,t) 
jf [G(x, z, t)E(z, y, t) - E(x, z, t)G(z, y, t)] (30) 



G^y,*) 



-G(x,y,t) 



+ 



2 / E(x,w,t)G(w,z,t)E(z,y,t) 



-i J [E(x, z, i)G(z, y, t) + G(x, z, f)£(z, y, *)] 



(31) 



because of We note that the density gets contributions from the condensate field ip 

as well as from the fluctuations (G,S). The contribution from ip*ip is the condensate density. 
We introduce the generalized potentials 



W d (x,y,t) = 5(x-y) jfp(z,z,t)V(x-z) 
U e (x,y,t) = p{^y ) t)V{-K-y)=Ul + iUl 
W p (x,y,t) = A(x,y,t)V(x-y)=Z4+iZ4 



(32) 
(33) 
(34) 



where the notation emphasizes real and imaginary parts of U p . We also define the matrices 



C(x,y,t) = /i(x,y)-p + Z4(x,y,t)+W d (x,y,*) (35) 

A(x,y,t) = C r (x,y,t)+^(x,y,t) (36) 

fl(x,y,i) = C r (x,y,0-^(x,y,0 (37) 
From Eq. (p6"l)-(|35|) we obtain an abstract matrix form of the equations of motion 



^G-'AG- 1 - - | + {Ul £} - [Ul S] 



(38) 
(39) 
(40) 



S = 

= aj>-Upij)*, 

where we have used the fact that (S,G) are symmetric matrices. These equations (p8|)-([40D 
are the nonlinear field equations for an arbitrary interaction V between the particles and 
contain any external potential through h. As an example the matrix product G~ l AG~ l can 
be written in coordinate representation as 



/ G- 1 (x,z,t)A(z,w)G- 1 (w,y,t). 

J z,w 



(41) 



The static equations can be obtained by setting the canonical momenta to zero, that is 
£(x,t,t) = 7r(x,t) = 0, G(x,y,t) = 0(x,t) = 0. From (H), © and (H) we than have 



i?(x,y,0) = i?(x,y) 
D(x,y,0) = D(x,y) 
^(x,0) = VW = 



1 
2 
1 
2 
fx) 



C-^^y) 



g _1 ( x > y) 



+ G(x, y) - <5(x - y) 



- G(x, y) 



(42) 
(43) 
(44) 



So that for the static case we have to self consistently solve 

~ / G- 1 (x,z)A(z,w)G- 1 (w,y)- J B(x, y ) = 

4 Jz,w 



£(x, z)^(z) - 2^(x)V> 2 (z)\/(x - z) 
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0. 



(45) 
(46) 



using fl3"2"|)-(|35|). We note that the constraint G = 1/2 for the static solution leads to R = 
D = and A = B so that equation ( |46|) becomes the usual non linear equation for the 



single quantity ip obtained from a many body permanent for bosons |fL2|| . However the time 
dependent eqs. (|38D-(|40D are more general, because our trial Gaussian is actually a coherent 
state with an indefinite number of particles. 

3 Uniform Static Solution 



To solve equations ( fig) and ( |4q) in the uniform case we go to the momentum basis where G is 
diagonal simultaneously with A, B and C . There exists a straightforward generalization for 
the nonuniform case where A, B and C do not commute. In the uniform case 0(k) = 05(k). 
The static equations ( f4"5| ) and (|46| ) become 

l -G-\k)A{k) - B{k) = (47) 
~B(0) -X(f) 2 ] = 0. (48) 

Where we have defined A = ^(0), V indicates the Fourier transform of the potential and 
e(k) = The diagonal matrices A and B (P5|)-(P7D can be written as 

A(k) = e(k)-/i + ^(k)+W c (k)+W d (k) (49) 
B(k) = e(k)-/i-^(k)+W c (k)+W d (k). (50) 

In this representation the generalized potentials (^)-(Q) become [we use the notation : 

U d (k) = X [ p(k') (51) 
W c (k) = / p(k')V(k-k') (52) 
ZY p (k) = / A(k')V(k - k') (53) 
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The solutions for (|47j ) and ((48|) are 



G(k) 



A(k) 



2\B(k) 
or B(0) = \<f) 2 , 



(54) 
(55) 



where the second solution is symmetry breaking because (ip) = (p/y/2 ^ 0. For stability 
reasons we shall see (in section 4) that both A and B must be positive which means that 
using (^) there is no ^ solution for A < 0. Using (|54]) with ([42]) and ([43]) we can express 
-D and R as functions of A, B 

1 5(k) - A(k) 



4 <M(k)2?(k) 



D(k) 



= l|B(k) + A(k)_ 1 



2 2,M(k)5(k) 



(56) 
(57) 



From (f28|), (|2D, (H) and (0) we have 



2 / £) k /y(k' - k) - A0 2 



2(e(fc) - /i) + A0 2 + / A + V(k - k') fl(k'). 



A(k)-B(k) 
A(k) + B(k) 

For a given density p of the system we can write the constraint relation M jV 
gives us 



(58) 
(59) 

p which 



(60) 

(||) and (§g). Note that 



p = <P 2 + f R(k>). 

To find a solution in the case <fi = we have to solve eqs. (^ 
.R and .D are functions of A, B through (pq ) and (|57|) . With these three equations one can 
obtain and /x. 

For the symmetry breaking solution <fi ^ 0, p can be determined using the equation 
B(0) = X(j) 2 which gives us 



/i = W e (O)+W d (O)-^(O)-A0 2 

Which means that solving (|58|)-(|6T|) gives us A, B and (f). 
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(61) 



4 Equations of Motion in the Small Oscillation Regime 
Generalized RPA 



To find the RPA equations we expand all the quantities around their equilibrium values, 
thus 2 

G(k,k',t) = G(k)5{k-k') + 5G{k,k',t) (62) 
E(k,k',t) = 5E(k,k',t) (63) 
4>(k,t) = <f>6(k) + 54>(k,t) (64) 
7r(M) = 5tt(M)- (65) 

For simplicity we write G and E in the basis where the equilibrium G is diagonal and keep 
terms up to first order in small quantities. For the uniform case the diagonal basis will be 
plane waves. It will be useful to introduce new momentum coordinates so that 

P = k-k' (66) 
k + k 

q = — ^— (67) 

so that 

5G{k, k') -> 5G(P, q) (68) 

and 

<fG*(P,q) = «?(-P,q). (69) 

Note that from flT8] ) the canonical variable to 5G(P,q,t) is 5E(— P,q). We will see that P 
and q can be interpreted as total and relative momenta respectively so that we can write the 
RPA equations in the form where P is diagonal and can be considered as a dummy variable 

8G(q, P, t) = s K (q, P)5E(q, P, t) + c K (q, P)5tt(P, t) + f S K (q, q', P)(JE(q', P, t)(70) 



'In equations ( |62[ ) and (|64|) the first delta is a Dirac delta and the second means an infinitesimal change. 
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-5±(q,P,t) = s M (q,P)5G(q,P,t) + c M (q,P)5<f>(P,t) + f 5 M (q, q', P)5G(q', P, $71) 

Jq' 

<ty(P,t) =57r(P,t)A(P)+ / c K (q , ,P)5S(q , ,P,t) (72) 
-57r(P,0 =8<f)(P,t)B(P)+ [ c M (q',P)5G(q',P,t). (73) 

>/q' 

Where we note that the (ir, 0) degrees of freedom are coupled to the much more numerous 
degrees of freedom (E, G) which are labeled by q. Different q values among (E, G) are 
also coupled. Introducing the notation /(q' + P/2) = and /(q — P/2) = /_, we find 
non-diagonal matrices in (q, q') 



S K (q,q',P) 
5 M (q,q',P) 



2F(q-q') [G + G" + + G_G"_ 
V(q-q') 



and diagonal elements 



+ 



2 

G^GZ 1 



+ 



GTG 



l/^-i 



t/(P) 



G' G\ 



V(q-q') 



(74) 



(75) 



s K (q,P) = 2[A h G_ + A_G + ] 
* M (q,P) = g 

Finally we see the coupling elements between (ir, <fi) and (E, G) 



c K (q,P) 

c M (q,P) 



G+ + G_ 



v+ + vc. + v"(P) i - 



(76) 
(77) 



(78) 
(79) 



which vanish when the symmetry is conserved (0 = 0). As pointed out above the equations 
are diagonal in P so we can interpret it as the total momentum of a pair of quasi-bosons, 
where — k' is thus interpreted as the momentum of a" hole" . Because 5G, 5E and 5(f), Sir are 
canonical variables we may invert the definitions of momentum and coordinate, for reasons 
that will be clear in the next section, where the zero mode is discussed. We define column 
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vectors 

(80) 

We can write a coupled oscillator Hamiltonian that corresponds to the RPA equations of 
motion in a matrix element form 

H RPA = ^(Q\A\Q) + \(U\B\U), (81) 

where the matrixes A and B are 

A -{ c K A J \ c M B J' (82) 

and the equations of motion are 

|0) = B|n) (83) 

|ri) = _^|e) (84) 

which can be written as second order equations 

|0) = -BA\e) (85) 
|fl) = -AB\U). (86) 

We may separate the diagonal part of H RPA so that 

Hrpa = Hq + H int , (87) 

where 

^>(o !)(£K<' G ' °b){%) ^ 

and ifi nt is equal to 
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Note that A and B must be positive for stability. If we introduce the trivial multiplicative 
canonical transformation 



5n) ~* \5ir) ~ \ Suy/B t 

5G\ (m: ' m; 



\S(f) 



&<t> 



(90) 
(91) 



we obtain 



Where ui and f2 2 are 



w(P) = V^( p ) 5 ( p ) ( 93 ) 



n 2 (q,P) = V s K(q,P)«M(q,P). (94) 

If we use the definitions of sk and sm from (|76|) and (|77"D we get after same algebra the 
remarkable result 



f2 2 (q, P) = \JA + B + + ^A.B_ = w(k) + u;(k'), (95) 

so that cu(P) and f2 2 (cu P) can be interpreted as the one and two free quasi-boson energies. 
We note that f2 2 (0,P) = 2u(P/2), which means that at zero relative momentum f2 2 (P ; 0) 
corresponds to two quasi-bosons each with momentum P/2. For the 7^ case the free 
quasi-boson energies u;(P) in general have a gap which means, that u(0) 7^ 0. In fact if we 
use and (§1]) we get 

A(0) = 2U r p {0) + \<f (96) 
B(0) = X(f) 2 (97) 

because of ( |29"| ) we write 

14(0) = J v A(k')V(k') = -\tf + I D(k')V(k'). (98) 
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We have then 



w(0) = yfA(0)B(0) = y 2A0 2 j£ V(k')£>(k') (99) 

The oscillations of the 50, 57r pair can be interpreted as a quasi-boson mode with a gap 
while the oscillations of SG, SH can be interpreted as an interacting pair of these same quasi- 
bosons. When = 0, we get S 3 = S e = and the one and two quasi-bosons systems are 
calculated independently. When ^ we must rediagonalize so that our final modes will be 
mixtures of one and two quasi-bosons. The variable q represents the internal motion of the 
quasi-boson pair with interaction given by the quantities S. In general this is a scattering 
problem and we must search for the scattering amplitude at a given energy and P, where 
the boundary conditions are determined by fl94|) . We can write the equations for SG and Scj) 
using 

SG(q, P,t) + / 5 K (q, q, P)«S M (q', q", P)5G(q", P, t) + c K (q, P) / c M (q', P)5G(q', P, t) 



5(P)c K (q,P)+ / 5 K (q,q / ,P)c M (q / ,P) 
Jo' 



6<f>(P,t) 



P,t) + 



A(P)B(P) 



CK(q , ,P)c M (q , ,P) 



P,t) 



(100) 
(101) 



-A(P) { c M (q\P)5G(q',P,t)- / CK (q / ,P)«S M (q / ,q // ,P)5G ! (q // ,P,t) 

Jo! Jq',q" 

where S contain both S and s defined in (|74] ) - (|77D i.e. 



<S K (q,q',P) = S K (q.(q'P) + ^(q ) q'P)%'-q) 
5 M (q,q',P) = S , M (q,q / ,P) + s M (q,q / ,P)5(q / -q). 



(102) 
(103) 



If we look for oscillatory solutions for the equations (100) and (101) 



5G(q,P,t) 
5E(q,P,t) 

Sn(P,t) 

(p,0 



SG(q,P,t = 0)e int 
5£(q,P,t = 0)e int 
7r(P,t = 0)e int 

<f)(p,t = oy nt 



(104) 
(105) 
(106) 
(107) 



14 



we may rewrite ( |100[ ) and ( |101| ) with the compact mode notation of (|86] 



(Q 2 - h 2 )U G = (c K B + J <S K ) n 
(Q 2 - h)^ = (aJcm + J c K S M ) n 



G ■ 



(108) 
(109) 



We have defined hi and h 2 as 



h x = AB + / c k cm 
h 2 = Sk<Sm + crCm 



(110) 
(111) 



where hi depends only on P and h 2 is understood as a P dependent operator in q and q' 
with ckCm being of separable form. We may write hi and h 2 in terms of the free quasi-boson 
frequencies defined in (|93|) and (p4|). 

(112) 
(113) 



hi = u" + I c k cm 
'q 



Q 2 + S'k'S'm + s k5m + SkSm = Q 2 + Vq. 



We note that, in this form, h 2 , ( 108|) and (|109|) are not hermitian. However it can be shown 
that they correspond to a fully hermitian form which leads to a unitary S matrix. We do 
not enter into this discussion here. | 

Eliminating IT^ from ( |108[ ) and ( |109| ) we get 

1 



(n 2 - hi) - (c k B + S K c M ) 



(Q 2 -hi+ ie) 



A J c M + J c K S M 



5G = 



(114) 



Defining the separable potential 

Vi = (c k B + S k cm) 



(VI 2 -hi + ie] 



(Ac M + c K S 



Ml 



(115) 



we can write 



(tt 2 -tt 2 2 ) = (Vo + Vi) 6G (116) 

3 If this were a problem with a discrete spectrum the non hermitian structure has the form AB with A 
and B both hermitian but non commuting. It is well known that the secular determinant Det|f2 2 — AB\ 
is equivalent to that of the hermitian form Det|f2 2 
variables II. 



AB"\A4| , corresponding to a transformation of the 
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which lead us to the standard form for the scattering for a given value for Q = ^(P, q) and 
^ = fi 2 (P,q'). 

II G (q, q') = tf(q - q') + ^— ^—VTl G = 5(q - q') + Q^VU G . (117) 

In the usual fashion we can define a T matrix as 

T = VU G (118) 

so that the equation that determines the matrix T is 

T = V + Vg { +) T. (119) 

In principle for a given interaction we can solve for the T matrix obtaining all the scat- 
tering properties like phase shifts and bound states. The existence of bound states will be 
determined by the poles of the T matrix. Because of the coupling between one and two 
quasi-bosons there can be one "bound" state which can be interpreted as the dispersion 
relation of a new RPA-boson Q^. Again using (|108|) and ( |109| ) we can eliminate SG in favor 



of 8(f) and obtain a formal form for the dispersion relation for f2t,(P) 

(til - hi) - (Ac M + c K S M ) Q2 l _ h (Bc K + c M S K ) = 0. (120) 

This involves the Greens function for h% below the threshold so that the ie in not necessary. 
Because of the broken symmetry (see next section) this dispersion relation will start at 
zero, in contrast to u;(P) given in (|9"5|). This can be seen schematically in Fig. 1. If the 
interaction is attractive enough there can of course be additional bound states. If any bound 
state actually leads to a negative value for Q 2 the system is unstable, as is usual for the 
RPA. 



16 



5 The Goldstone mode 



In order to understand the structure of the Goldstone mode we recognize from d85|) that for 
Q = we must have 

A\e m ) = o, (121) 

going back to the more explicit notation we have 

/ 5K(q,q / ,P)4 0) (q / ,P) + ^(q,P)4° ) (q,P) + CK(q,P)^ 0) (P) = (122) 
/ c K (q' J P)0j? ) (q , ,P) +4 0) (P)^(P) = 0. (123) 

J a' 



Eliminating #[ 0) (P) from ( |T2"BD we get an equation for 6 



(0) 



£ 5 K (q, q', P)4 0) (q, P) + *K(q, P)4 0) (q, P) - ^7pJ<*(q, P) jf <*(q', P)4 0) (q', P) = 0. 

(124) 

If P = we can find a solution for (|124 ) 







(0) 



(125) 



which arises from the properties of the static solutions. To do so we substitute ([125D in fll2l[) 
and using the static properties fl5"8|)-(|5T3D we compute each one of the integrals separately. 
For convenience we omit the index (P = 0) 



q' 



%(q, q') 



G- 2 (q) 



4G(q) / v(q-q) 



G(q') 



4G(q) S(q) - A(q) - <f> 2 V{q) 



(126) 



and 



CK(q 



G' 2 (q) 



20 / V(q') 



G(q') 



20 5(0) - A(0) - A0 2 = -20,4(0). 



(127) 
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It is easy to see that the subtitution of ( |126| ) and ( |1 27| ) leads to satisfying ( |124| ). Once we 
know 9jP we can find 9^ using (|123|) so that 



9^ 



A(0) 7 q 



CKlq 



^ 2 (q) 



128) 



This gives us a particular normalized column vector 9^ for the zero mode 



= c h 



G- 2 (q,0) 
4 



with 



C 



4' \ 



G- 2 (q,0) 



(129) 



(130) 



+ 



This corresponds to a new coordinate and momentum which are linear combinations of the 
previous (SG, S(p) and (<5E, Sir) and can be written as inner products. 



Q = (0<°>|e) = c(± 

V = (^°)|II)=C(| 



1 - 



G- 2 (q,0) 



G" 2 (q,Q) 



5G 



(131) 
(132) 



To better understand these coordinates we expand Af, the mean number of particles, up to 
first order 

Af = Af + 5N, (133) 



where A/q is the static result given by 



+ G(q',0)-1 



(134) 



The first order term 8J\f is 



Which we see is proportional to our new canonical momentum 

V = C5M. 
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(135) 



(136) 



From (HD and ( |P2T|) we have 

(0(°) |li) = P = -(# (0) |^|0) = 0. (137) 

Indeed this is the general result (5Af = 0) proved in section 1 and, as expected, the zero 
mode corresponds to a "translational" motion in N with no oscillation. 



In fact because we have the RPA Hamiltonian ([y]) it is possible to see the explicit form 
for the mass coefficient M. in the quadratic expression for the 8M dependence. To see this 
we first rewrite the RPA Hamiltonian ( [Sip in a basis where A is diagonal, because it contains 
the zero mode. 

#rpa = \f{i\A\x)# + \p v °(i\B\j)P : ''■ (138) 

with 

(0|.A|0) = 0. (139) 

This shows that there is no dependence on 9^ which corresponds to our original parameter 
9 canonical to TV. The summation on is implicit. We now separate the contributions 
p(o) from the kinetic piece as 

K = ~P** {i\B\j)P j = ^P m *(m\B\m')P m ' + ^P°\0\B\m)P m 

+ip m * (m\B\0)P° + \p°* (0|B|0)P° (140) 

where m and vnl are different from zero. We can diagonalize B in the subspace getting 

K = ip A *(A|£|A)P A + ip°>|£|A)P A + ip A *(A|£|0)P° 

Zi JU JU 

+^P°*(0|i3|0)P°. (141) 
Using the fact that P° is constant we complete squares 



p \ + p o x x 



(\\B\X) [P x + X X P°] + ^P°*(0|£|0)P° 



-ip°*X A *(A|£|A)X A *P° (142) 
2 
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where the vector X has to be equal to 



• YA = (0|e|A) (A^A)- (143 > 



So that the coefficient of (P°) 2 /2 will be 



= (0|B|0) - (0|B|A) p^y(A|^|0) (144) 



which gives us 



M (0\B\0) - (0|S|A)^(A|S|0) " ( % |0) + (0|B|0)- (M5) 

Note that this a non-trivial result, and it is the effect of the linear terms in P° in (|140|) . If 
we use ( |136| ) we can actually see that 

2M 2M K J 

which means we have calculated the coefficient of 5M 2 . To see the physical meaning of this 
coefficient we expand our original Ti — fjj\f around the equilibrium value JVq 



h{N) - t*VW = + SAf-^r + \^ 2 ^i - + ■■■■ (147) 



Using 



we have 



" = M < 148 > 



H(N) - »{N )N = H(K) - + \ M2 ^m ( 149 ) 



comparing the 5N" 2 coefficients in ( |146| ) and (|149 ) we have 



d 2 H dfj, C 2 



(150) 



dAf$ dAf M 

Which means that once we have computed KA using ( |145|) we also get the value of the second 
derivative of the energy with respect to the number of particles. 
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For the remaining modes we have to work with a Hamiltonian in the sub-space that 
excludes the zero mode i.e. 



(151) 



If we introduce the canonical transformation 



(152) 
(153) 



we get 



#rpa = \pfPl + \<£ [ <\ \IB*A»\IB»\.i ) Q<- 



(154) 



If we now diagonalize the matrix y/B^A s y/B^ we have the final form to be used in subsection 
6.2. 

H S RPA = lp n *V n + ^Q n *n 2 Q n , (155) 



where Q 2 is real for stability and the V n and Q n are in general complex. 



6 General Remarks 

6.1 Connection with the Bogoliubov Transformation 

One can construct the operator that annihilates our trial wave functional (H) so that 



£(W,t] = 



(156) 



where 



and 



{(*) = 5 



Ql/2 



4iG 1/2 T, 



(157) 



571 = 7T — 7T. 



158) 
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Note that this result is valid for any particular time t. For the present purpose we will use 
the static annihilation operator 



1(0) = i = + iG 1 ' 2 



71 



(159) 



2GV2 

where now both G and <fi are the time independent equilibrium results. This will be our 
quasi-boson annihilation operator. The corresponding creation operator is then 

1 1 



and as usual 



2GV2 



iC 1 ' 2 * 



(160) 



a.b 



K.,eb] = [^,S] = o 



then 



and 



7T 



2VG 



(161) 
(162) 

(163) 
(164) 



We can now make the connection with the Bogoliubov transformation writing the quasi- 
boson operator function of the original ip as 



Z = x ip-ip 



+ Y 



(165) 



with 



X 
Y 



1 

2^2 



Ql/2 
1 



+ 2G X ' 2 
- 2G 1 ' 2 



(166) 
(167) 



Equation ( |165| ) is a generalized version of the Bogoliubov transformation where the commu- 
tation rules (|161| ) imply that the matrix relation 



X 2 - Y 2 = 1 
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(168) 



as 



D = XY (169) 
R = Y 2 (170) 

We may write the original Hamiltonian in terms of these quasi-boson operators £ and £t. 
It is easy to see that there will be terms like (friH^ which turn two quasi-bosons into one 
justifying our interpretation of the generalized RPA. In the next subsection we will show how 
the small oscillations correspond to modes that are mixtures of one and two quasi-bosons. 

6.2 One and two quasi-boson coupling 

We proceed by expanding our wave functional around the equilibrium values for G, E, <p and 

7T 

(0'|*(t)|> = 9[<f/,t] « [l - J 5<t>' (~5G± - *<5s) 5<P' - J 5<p' [^5<p - z5tt)] * [^ 

(171) 

where [</>'] is the ground state wave functional that is annihilated by £. If we use 



and (|164j) we can write 

w, t] « |i - / \~=sg^= - iVGssVG^ + a + ee f + 

^=5<p - iVGS-K^j + } *oM- (172) 



Note that the equilibrium G can be taken diagonal in a particular representation while in 
this same representation SG and <5E are non-diagonal. If we expand 5G , 5(f), <5E and for in 
the Q modes found by solving the RPA equations we have the general structure. Thus 

5E\ = ^(%)^ (173) 



for 
5G 



(*g\ ( 174 ) 
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with the hamiltonian ( 155|) for V^ 1 and Q n . The state vector looks like 

l*W>«{|0>+£[|l> + |2>]e«*} (175) 
where we have taken into account that £ annihilates the ground state and have defined 



|1> = J*??\0) (176) 
|2) = [Q$M\0) (177) 



where $^ and $2 are the one and two particle wave functions given by 

$ n = <p n -^=6^ -iQ n ^G0^ (178) 
2v G 

$ n = pn_*n_^_ i gn v ^0n > /£ i ( 179 ) 
2vG 2vG 

Note that for each value of Q (either the bound state uj or the the continuum Q 2 ) we have, 
as expected, a mixture of one and two quasi-bosons. We take into account that for each 
given mode [] have a relationship between V and Q n 

v n 

Q n = (180) 
Then the unnormalized wave functions can be written as 



*? = ^m%-^§-% ( 181 ) 



$2 = — W^— U - ^VGe^VG. (182) 

In order to understand the significance of the bound state for P = we can examine its time 
independent two body component, $^(P = 0,q). In this particular case 

$° = (183) 



4 In fact as we note before in section 5 this relation does not hold for the zero mode where V = constant 
and Q° = (V/M)t 
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From ( |125| ) we get 

$9 = _ 

If it is to represent the bound state of a pair of quasi-bosons the integral 



*2 = "4= ( 1 - TT^TT ) -4=^- (184) 



q 



B(q) 



^(q) - g(q) 



(185) 



/ q A(q) 

must be finite. To show this, we examine the behaviour of ($2) 2 . For q — > 0, it goes to a 
finite value because of the gap as discussed in fl99|). For q — > oo all the generalized potentials 
go to zero and both A and B go to infinity as q 2 plus a constant. This leads to {&%) 2 — > a/fc 
where a is the difference of these constants. This behavior assures the convergence of ( |185| ). 

For the case when P / we expect this normalizable bound state to evolve with P and 
to exist for at least a finite range of P. We note that it is not separable in P and q so that 
the relative wave function changes with P and could become unbound for some critical P. 

7 Temperature Dependent Calculation 



Temperature is often introduced [13[] in the so called temperature dependent Hartree Fock 
Bogoliubov approximation by generalizing the static quantities R and D defined in (f42|) and 
(||) as 

D(P) = X(P)Y(P) -> X(P)F(P)(1 + 2i/(P)) (186) 
R(P) = F 2 (P) ^X 2 (P)z/(P) + F 2 (P)(l + z/(P)) (187) 

with the occupation z/(P) given by 

^( p ) = e , ( P)/L _ i ( i8s ) 

This corresponds to the usual grand canonical ensemble where u is the one quasi-boson 
energy defined in (|93l) , with its gap. This obviously describes the physics as having free 
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quasi-bosons where neither the existence of the gapless bound state nor the interaction 
between the quasi-bosons is taken into account. 

Our approach to the problem is different in that we use the canonical ensemble. We have 
seen that the RPA Hamiltonian consists of a free quasi-boson part H Q ( |8~8"D and H int (|89|) that 
treats the interaction between the quasi-bosons. So our first approximation in the canonical 
ensemble for the free energy will be 

F«(A/\ T)=U + \I log(l - e-^ p >). (189) 
p Jp 

where u is the free quasi-boson energy given in fl93|). We can also extract more information 
due to the interaction of the quasi-bosons and this can be done by computing the second 



virial coefficient 10 



F (2 \fV,T) = -J e-^( p ) + — / e-^TrS-'^-S. (190) 
p Jp Airi Jp, q oih 

Where O(P) is the zero mode bound state dispersion relation and S can be determined, 
as mention before, through the scattering of the quasi-bosons using ( |119| ). Note that for 



the cases where we have two variational solutions = and <p ^ we can calculate using 
( |189| ) and ( |190| ) for a given M obtaining two free energies and this allow us to determine the 



critical temperature by setting their derivatives with respect to Af equal i.e looking for the 
same value of the chemical potential [O] 
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Figure Caption 

Fig. 1: Schematic plot of the one u and two ^2 free quasi-boson as a 
function of the total momentum. When the interaction between them is 
taken in account (RPA) we have a scattering region and a gapless bound 
state Q 
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